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Abstract 



As is well known, structure formation in the Universe at times after de- 
coupling can be described by hydrodynamic equations. These are shown here 
to be equivalent to a generalization of the stochastic Kardar-Parisi-Zhang 
equation with time-dependent viscosity in epochs of dissipation. As a conse- 
quence of the Dynamical Critical Scaling induced by noise and fluctuations, 
these equations describe the fractal behavior (with a scale dependent fractal 
dimension) observed at the smaller scales for the galaxy-to-galaxy correla- 
tion function and also the Harrison-ZeFdovich spectrum at decoupling. By 
a Renormalization Group calculation of the two-point correlation function 
between galaxies in the presence of (i) the expansion of the Universe and (ii) 
non-equilibrium, we can account, from first principles, for the main features 
of the observed shape of the power spectrum. 

Subject headings: Gravitation; Cosmology; Critical Phenomena. 
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Contemporary cosmology is riddled with problems like the Problems of the Cosmological 
Constant, Dark Matter and the Age of the Universe. In addition, we have no dynamical 
understanding of why the two-point correlation function for galaxies, a key point of contact 
between theory and observation, has the observed behavior. It is known from observations 
that the galaxy-to-galaxy correlation function, £oss( r )> * s wen fit by a power law of the 
form 

£oBs(r) oc r~ 7 , 

where r is the comoving separation between the galaxies and 7 is determined from catalogs 
to be between 1.5 and 1.8 at distances of the order of the megaparsec. Nevertheless, at 
the epoch of decoupling, we know from COBE data Q that the correlation function goes 
like r -( 4 2±0 ' 3 ). How is this so? How can this deviation from an integer in the power law 
exponent be accounted for? Why is there an evolution in the exponent of the power law? 
Can this fact be established from some generic physics? 

In this letter we will compute the galaxy-to-galaxy correlation function from the hy- 
drodynamics that describes the formation of structure in the Universe and will present an 
answer to the above questions. 

We generalize previous results obtained by Berera and Fang in Ref. and, independently 
by the authors of Ref. ||; we also compute the power spectrum of density perturbations. 
In the first case we have generalized their calculation to include the effects of self-gravity, 
expansion and, simultaneously, non-equilibrium; we have also generalized their asymptotic 
calculation to the full range of distances, from COBE and well into the realm of the galaxies. 
In the latter case, we have included deviations from equilibrium and the expansion of the 
Universe also in a self-consistent way, and we have extended their calculation back into the 
decoupling era. We will not need to introduce any new physics: our conclusions follow solely 
from a straightforward (albeit non-naive) analysis of the hydrodynamic equations and the 
extension to the Dynamical Renormalization Group and Dynamical Critical Phenomena of 
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techniques familiar from Condensed Matter and Elementary Particle Physics. 



As argued in Refs. and JJ, if the power law behavior of the two point correlation 



function for galaxies is due to some form of critical phenomena, it follows that in the realm 
of the galaxies there must exist some kind of fluctuations which should account for the 
observed behavior of £(r). We model them by means of power law correlated noise. 

The suitably averaged value of the 2-point correlation function for the density contrast, 
8(f,t),, written in comoving coordinates, is identified in phenomenology with £,oBs{ r ) (see, 
e. g., references (§], || and |§). We will study the scaling behavior of the contrast-contrast 
correlation function. 

Under the assumption^ of irrotational peculiar velocity u and peculiar acceleration w, it 
is straightforward to check that the hydrodynamic equations arising from the application of 
Newtonian considerations to structure formation in the Universe, can be written in comoving 
coordinates in terms of a velocity potential ip and the gravitational potential due to the 



contrast as (Ref. §,§,§) 



Of 



■V + Hi/, - ^(Wf --</>+ -f[p b {l + (4vrGa 2 p fe )- 1 V 2 0}] = 



(1) 



— (j) + H<j> - AnGap b ilj + T = 



(2) 




(3) 



u = —Vtjj 



(4) 



w = —a V0 



(5) 



V 2 = AitGa 2 Pb d(r, t) 



(6) 



x The rotational components of the velocity decouple at a quicker rate than their non-rotational 
counterparts, and thus for late enough times it is always possible to justify this assumption. 
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p(f,t) = p b [l + 5(r,t)}. 



(7) 



Here, a(t) is the scale parameter of the homogeneous cosmological background, and the 
function J-(r, t) originates in the gauge freedom associated with the irrotational characters 
of u and w. H is the Hubble parameter. The function / (also known as minus the specific 
enthalpy) is determined by the equation of state assumed for the matter whose clustering is 
described by the above equations: p = f(p), f'(x) = —f'(x)/x. We know however that for 
dust (p = 0), the Zel'dovich approximation, w = F(t)u, with 



works very well during the early non-linear regime. One can (and we will) assume that this 
approximation is also valid when pressure is included and provided that the pressure is small 



F = AiiGpb — H F — F 2 



(8) 



(cf. Ref. i). 



Under the above assumptions we have 




(9) 



S(x, t) 



lD(t) 



V 2 ^(f , t) 



(10) 



ai)(t) 



where D(t) is the growing mode component of the density perturbation (Refs. [ffl, |]5| and 
0). The velocity potential satisfies 



- ^m? = ^G Pb ^ + i c >g fi + -S^vV) , (ii) 

at a 2a D(t) a \ a D(t) J 



and expanding the logarithm to lowest order, we get an equation of the form 



0_ 

dt 



(12) 



where the fi(t) are determined by the background geometry as 



fi(t) 



c 2 s D(t) 
a 2 (t)D(t) 



(13) 
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m = UT) < 14) 
/s(() . toGp6(() ||_|| (15) 

This is a generalization with time-dependent coefficients (and a mass-like term) of the 
Kardar-Parisi-Zhang (KPZ) equation 

JU = y^ h + \\{Whf + T](X, t) , (16) 

which plays a central role in surface growth phenomena, and whose scaling behavior in the 
IR (large distance, small k) and UV (short distance, large k) regimes are well understood 
in terms of the correlation properties of the noise or in the absence of noise, the correlation 
properties of the initial conditions (e.g., Ref. ||). 

By a series of changes of variables and rescalings, one can rewrite equation (|T2|) as 



^-H(x, r) = V 2 H(x, t) + (Vff(x, r)) 2 + 9A ^ 9t H(x, r) + r,(x, r) (17) 

with A{t) = ^|^|exp ^I dr' and the noise term is rj = (/ 2 / fi)fj- The quantity H(x,t) 
is related to ip(x,t) by H(x,t = J T dr' fi(r')) = j^tp(x,t). The first term in the rhs 
represents the smoothing effect of diffusion and the second term is due to non-equilibrium 
effects. The third term describes the effects due to the time-dependence and matter content 
of the background geometry (it contains effects from the expansion of the Universe and 
self-gravity); it behaves, in the linear approximation and in a variant of "conformal-time" , 
as a mass term, and therefore introduces a natural correlation length into the problem. 
Finally, the noise term models the various fluctuations that can appear during epochs in the 
evolution of large scale structure in the Universe. 

It is clearly seen from Eq. ([T7| ) that the scaling properties of H(x, r) depend on two key 
features of the equation: (a) the characteristics of the noise and/or (b) the specific features 
of the background geometry. 
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The power spectrum P(k; t) is the Fourier transform of the two-point correlation function 
for the density contrast (see ||], and ||). We get0 

P(k;t) = <S>(t)k 2 q 2 (H(k,r) H{q,T)) dassical \- = _ % (18) 

where the two-point function for H(k, r) does not include the effects of non-linearities nor 
higher order effects due to fluctuations, i. e., it is computed from the classical theory, and 
$ 1 / 2 (t) = dD{t)/dt comes from the changes of variables needed to transform from Eq. 
(1) to Eq. (0). 

The renormalized two-point correlation function (H(k,r)H(q,r')) obeys a Callan- 
Symanzik equation, Ref. 0, the solution of which contains all its scale dependence, and 
P(k; t) can be written as 

P(k; t) Full = <$>(t)k 2 q 2 (H(k, r) H(q, r)) Improved (19) 

where (H(k, r)H(q, t)) Improved is the solution to the Callan-Symanzik equation. As is done 
within the context of the Renormalization Group (RG), this object is computed by inserting 
into the free two-point correlation functionQfor (H(k, T)H(q, r)) the values of the couplings 
obtained by solving their RG equations. 

We will now restrict ourselves to the following cosmological and noise scenario: flat FRW 
cosmologies where the noise is arbitrarily power-law correlated both in space and in time 0. 
In fiat FRW it is straightforward to see that $ FRW (t) = (81/4)c^o /3 ^ 4/3 - 

For noise with properties given by 

( V (k,uj)) = (20) 



2 In what follows and to simplify the writing, we will write down H(x,t) instead of H(x,t 
fdr'hir')). 
3 Obtained from Eq. (|l6|) by setting A = 0. 
4 A cosmology with constant A is flat FRW. 



{ri(k,u)Ti(q,n)) = 2 



D + D e k~ 2p 



UJ 


-29' 


LO 





5 3 {k + q)5{uj + n) 



(21) 



a straightforward calculation gives (here for convergence of some integrals, — 1/2 < 9 < 1/2) 



(H(k,r) H(q,r)) Improved 



v 2 (k) f U {k) U e {k) sec(irl 



S 3 (k + q) 



(22) 



27rir 3 A 2 (A;) [ k 2 ' k 2i - 1+2e +p) 
where we have introduced U = ^-K 3 , Ug = ^-K 3 , and K 3 = 

The running coupling constants v(k), X(k), Uo(k) and Ug(k) obey the following renor 
malization group equations (Ref. jlPfl), 

du 



-/i— = v 
dp 



~T2 Uo + 2 ~ R yT u ^ 1 + 26) sec(7r6 



^ = Ia^(1 + 29) sec(7r6 
a/i 3 



(23) 



(24) 



dUo 
d\i 



{49 - 1 + 2p)C/ e + 



3 + 20 
12 



We 



3 + 100 - 6p - 4flp 
12 ' 



1 + 29) sec(n9)U$ (25) 



-U n + -U 



^ + 9 6p + 80 (1 + 2Q) sec ( n9)UoU() + 1^2 (1 + 40) sec(27r0) (26) 



The parameter yu has dimensions of momentum. These equations have several fixed points, 
and the asymptotic behavior of the solutions depends on the values of the noise parameters 
9 and p (see Fig. [I] for a plot and definitions of the parameter regions.) The fixed point 
structure in the region of the (Ug — Uq) plane where the correlations are positive is shown 
in Figure 1. 

^From COBE observations we know that in the IR (small fc-regime) the power spectrum 
is Harrison-Zel'dovich; similarly at "large" momentum, the power spectrum is also scale 
invariant but with a different exponent. Therefore the boundary condition on the improved 
two point function is that for A; — » the power spectrum be Harrison-Zel'dovich, i. e., 
limfc_>o PFuii(k'i t) ~ k. This translates into a condition that must be satisfied by the two 
noise exponents 



2 - 49 - 2p = 1 
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(27) 



the subsequent evolution of P(k) is controlled by the RG according to the above equations. 
The results of integrating the RG equations in ( p3| ) - ( p6|) with a typical set of initial 
conditions compatible with COBE data are shown in Fig. 0. They reproduce both the 
general form of the power spectrum and its main observational features. However, as one 
evolves into shorter distances (larger momenta) higher order terms in the expansion of 
the logarithm in the hydrodynamic equation must be included since they become relevant 
at these shorter scales; the behavior displayed by the power spectrum in the calculation 
presented here, indicates that the evolution when higher order terms are included will also 
be in the direction of the observational evidence, since the largest momenta at which our 
calculation can be trusted is also consistent with the behavior inferred from catalogs of 
galaxies. 

We have shown that the hydrodynamics of a fluid of galaxies interacting through gravity 
can be studied using scaling techniques based on the dynamical renormalization group. We 
have taken into consideration all the effects due to self-gravity, expansion of the Universe 
and non-equilibrium present in the "fluid of galaxies". We have applied these ideas to 
the calculation of the 2-point correlation function for the density contrast, and have found 
that (i) its scaling behavior depends on the background geometry and the noise and/or 
initial conditions for the density contrast, but (ii) can be computed and (hi) comparison 
of our results with observations shows excellent agreement. In summary, we have seen 
that the power spectrum can be viewed as evidence of dynamical critical behavior in the 
Universe. In fact, because of this critical behavior, once the power spectrum at decoupling 
is known to be Harrison-Zel'dovich what happens at smaller scales is fairly insensitive to 
small deviations from this initial condition, since criticallity implies that the system will 
eventually be attracted to one of its fixed points irrespective of the details of the physics. 

Although we have demonstrated the feasibility of our approach for the simplest back- 
ground (flat FRW), it is clear that more general and complete cases can be similarly treated. 
Furthermore, these considerations lead to a very interesting view of the large scale structure 
of the Universe, where all kinds of new phenomena and behaviors can now be described; these 
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include pattern formation, roughening transitions, nucleation, defect generation, ecological- 
like behavior for galactic and other many body gravitational systems, etc.. 
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FIG. 1. Noise spectrum and fixed points for Uq and Uq. The fixed point structure for these RG 
equations depends on the noise spectrum. There exist five distinct regions in noise parameter space 
leading to three un-related behaviors of the RG equations. In the upper panel we display the noise 
"space" and in the lower panels the corresponding map of Uq and Uq fixed points together with 
their IR or UV characters. Different choices in the p~9 plane lead to qualitatively different solutions 
to the RGEs for Uq and Uq. In the region where correlations are positive, different behaviors can 
be found as shown in this figure. The difference between I and V (or the set III, IV, VI and VII) 
appears in the region where Uq < and where the RGEs have different types of critical points. 
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FIG. 2. The predicted power spectrum for the contrast as a function of the momentum scale 
and some typical observational data. The data plotted in the figure are taken from the CfA-101 
and CfA-130 catalogues (Ref. [10]). The steep fall-off in our predicted curve past its maximum 
indicates the need to include higher order terms in the expansion of the logarithm in Equation (11) 
which, as explained in the text, become relevant at shorter distances. 
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